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Abstract. Boundary value problems for integrable nonlinear evolution PDEs for- 
mulated on the half-line can be analyzed by the unified method introduced by one 
of the authors and used extensively in the literature. The implementation of this 
general method to this particular class of problems yields the solution in terms 
of the unique solution of a matrix Riemann-Hilbert problem formulated in the 
complex fc-plane (the Fourier plane), which has a jump matrix with explicit {x,t)- 
dependence involving four scalar functions of k, called spectral functions. Two of 
these functions depend on the initial data, whereas the other two depend on all 
boundary values. The most difficult step of the new method is the characteriza- 
tion of the latter two spectral functions in terms of the given initial and boundary 
data, i.e. the elimination of the unknown boundary values. For certain boundary 
conditions, called linearizable, this can be achieved simply using algebraic manip- 
ulations. Here, we present an effective characterization of the spectral functions in 
terms of the given initial and boundary data for the general case of non-linearizable 
boundary conditions. This characterization is based on the analysis of the so-called 
global relation, on the analysis of the equations obtained from the global relation 
via certain transformations leaving the dispersion relation of the associated lin- 
earized PDF invariant, and on the computation of the large k asymptotics of the 
eigenfunctions defining the relevant spectral functions. 
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1. Introduction 

A unified method for analyzing boundary value problems, extending ideas of the 
so-cahed inverse scattering transform method, was introduced in [29 \ \30 \ 138], This 
method has been implemented to linear and integrable nonlinear evolutions PDEs 
on the half-line and the finite interval PElElIISlIigiSlEaEIlEaEllESlESl 
EiliailHlISOlIMlETlEIlEaiMliniEaE^ linear and integrable 

nonlinear evolution PDEs in two space variables [32l [331 133 EH EH] , to linear elliptic 
PDEs [21 El H [El EH ia Sg Ea Eg EI], and to the two prototypical integrable 
nonlinear two-dimensional elliptic PDEs, namely the sine-Gordon and Ernst equations 

[5I1E21E31E51E1]. 

Regarding linear evolution PDEs containing spatial derivatives of second order, the 
unified method yields novel integral representations which have both analytic and 
numerical advantages in comparison with the classical integral and series representa- 
tions: (i) They are uniformly convergent at the boundaries; analytically, this makes it 
easier to prove rigorously the validity of such representations without the a priori as- 
sumption of existence, whereas numerically, using appropriate contour deformations, 
it makes it possible to obtain integrands which decay exponentially as \k\ — )• oo and 
this leads to efficient numerical computations \27\ I70j. (ii) These representations re- 
tain their form even for more complicated boundary conditions, whereas the classical 
representations involve infinite series over a spectrum determined by a transcenden- 
tal equation. For example, in the case of the heat equation with Robin boundary 
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conditions, the classical representation involves an infinite series over {A;„}j^ where 
kn satisfy a transcendental equation, whereas the new method yields an explicit in- 
tegral representation. For linear evolution PDEs containing x-derivatives of arbitrary 
order formulated on the half-line or the finite interval, the unified method still yields 
integral representations explicitly defined in terms of the given initial and boundary 
conditions, whereas it is shown in [69j and [78j that even for the linearized KdV on 
the finite interval with generic boundary conditions, there does not exist an infinite 
series representation. 

Regarding integrable nonlinear evolution PDEs in one space variable, the unified 
method yields novel integral representations formulated in the complex fc-plane (the 
Fourier plane). These representations are similar to the integral representations for 
the linearized versions of these nonlinear PDEs, but also contain the entries of a 
certain matrix- valued function M{x,t,k), which is the solution of a matrix Riemann- 
Hilbert (RH) problem. The main advantage of the new method is the fact that 
this RH problem involves a jump matrix with explicit (x, t)-dependence, uniquely 
defined in terms of four scalar functions called spectral functions and denoted by 
{a{k), b{k), A{k), B{k)}. The functions a{k) and b{k) are defined in terms of the initial 
data ^0(2;) = q{x,0) via a system of linear Volterra integral equations. The functions 
A{k) and B{k) are also defined via a system of linear Volterra integral equations, but 
these integral equations involve all boundary values. For example, for the nonlinear 
Schrodinger (NLS) and the modified Korteweg-de Vries (mKdV) equations formu- 
lated on the half-line, A{k) and B{k) are defined in terms of {q{0,t),qx{0,t)} and 
{q{0,t), qx{0,t), qxx{0,t)} respectively. A major difficulty of initial-boundary value 
problems is that some of these boundary values are unknown. For example, for the 
Dirichlet problem of the NLS and of the mKdV, the boundary values qx{0,t) and 
{qx{0,t), qxx{0,t)} respectively, are unknown. It turns out that this difficulty can be 
bypassed by utilising the so-called global relation, which is a simple algebraic equation 
which couples the spectral functions. 

1.1. Linearizable versus non-linearizable BVPs. We now distinguish two cases: 
(a) For a particular class of boundary conditions called linearizable, it is possible to 
express A{k) and B{k) in terms of {a{k),b{k)} and the given boundary conditions. 
This can be achieved by analyzing the global relation and the equations obtained 
from the global relation under those transformations in the complex fc-plane which 
leave invariant the dispersion relation of the linearized version of the given nonlinear 
PDE. It must be emphasized that the above analysis is carried out in the A:-space, 
thus the spectral functions A(k) and B{k) are determined directly, without the need 
of determining the unknown boundary values. In summary, for this class of problems 
the unified method is as effective as the classical inverse scattering transform method. 
(6) For non-linearizable boundary conditions which decay for large t, by utilizing the 
crucial feature of the new method that it yields RH problems with explicit {x, t)- 
dependence, it is possible to obtain useful asymptotic information about the solution 
without the need to characterize the spectral functions A{k) and B{k) in terms of 
the given initial and boundary conditions. This can be achieved by employing the 
Deift-Zhou method [22l [23] for the long-time asymptotics gOl iH ill |l3] and the 
Deift-Zhou-Venakides method [2(Hl21j for the zero-dispersion limit |46[|59j . However, 
the complete solution of non-linearizable boundary value problems requires the char- 
acterization of {A(k), B(k)} in terms of the given initial and boundary conditions. 
The effective solution of this problem is particularly important for the physically sig- 
nificant case of boundary conditions which are periodic in t, since in this case it is not 
possible to obtain the rigorous form of the long time asymptotics, without the full 
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characterization of A{k) and B{k) (in spite of this difficulty important resuhs about 
the large t-asymptotics of such problems are presented in [lOl [HI [12] ) . 

In this paper, we revisit the implementation of the unified method to integrable 
evolution PDEs on the half-line. We concentrate on the effective characterization of 
the spectral functions for non-linearizable boundary conditions. We call a character- 
ization of {A[k), B{k)} effective if it fulfills the following requirements: (a) In the 
linear limit, it yields an effective solution of the linearized boundary value problem, 
i.e. it yields a solution in the form of an integral which involves the transforms of the 
given initial and boundary conditions. (6) For 'small' boundary conditions, it yields 
an effective perturbative scheme, i.e. it yields an expression in which each term can 
be computed uniquely in a well-defined recursive scheme. 

The effective characterization of {A{k)^ B{k)} presented here is based on an effec- 
tive characterization of the unknown boundary values in terms of the given initial and 
boundary conditions, i.e. on an effective characterization of the so-called generalized 
Dirichlet to Neumann map. An effective characterization of this map for the Dirichlet 
problem of the NLS was presented in [7j by employing the so-called Gelfand-Levitan- 
Marchenko (GLM) representations. Following this important development, analogous 
results for the Dirichlet problem of the sine-Gordon, mKdV, and KdV equations were 
presented in [35', 'HS] . The final formulas of [35] did not involve the GLM representa- 
tions, but their derivation was based on these representations. Here, we rederive the 
formulas of |35j for the Dirichlet problem of the NLS and of the mKdV directly, i.e. 
without using the GLM representation, and also present the analogous formulas for 
the Neumann problem of the NLS as well as for the mKdV equation in the case that 
either (/x(0, t) or 3^re prescribed as boundary conditions. 

Our approach uses three ingredients: (a) The large k asymptotics of the eigenfunc- 
tion <I>(t, k) defining {A{k)^ B{k)}. (b) The global relation and the equations obtained 
from the global relation under the transformations which leave invariant the dispersion 
relation of the associated linearized equation, (c) A perturbative scheme to estab- 
lish effectiveness. We emphasize that, in general, all three ingredients are needed, 
except for the particular case of the NLS, where as explained in appendix [A] it is not 
necessary to use the invariance of the global relation. In particular, ingredient (a) 
yields several possible formulas for the unknown boundary values, so it is absolutely 
necessary to employ ingredient (c) in order to choose the one that yields an effective 
solution. 

The important idea that the asymptotics of the associated eigenfunctions can be 
used for the characterization of the unknown boundary values for the NLS was first 
introduced in [T8] (the authors of [18j employed the asymptotics of the eigenfunction 
F{t,k) instead of ^{t,k), but these eigenfunctions are related by equation (2.16)). 
In other words, the ingredient (a) above was first introduced in [18] . However, the 
ingredients (b) and (c) have not been introduced before. Hence, the approach of 
|18] cannot be generalized to other integral nonlinear PDEs, such as the mKdV. 
Furthermore, since the ingredient (c) was not introduced in [18j . it is not a priori 
clear if any of the formulas presented in [18] yield an effective characterization. It 
turns out, as explained in appendix |XJ that the relevant formulas for the NLS are 
effective in the case of the half-line but not for the case of the finite interval ^67j . 

From the above discussion it follows that it is possible to construct an effective char- 
acterization of both the generalized Dirichlet to Neumann map and of {A{k), B{k)}, 
without the need to introduce the GLM representations. However, it appears that the 
latter representations might have an advantage for analyzing the large t-asymptotics 
of the solution in the case of t-periodic boundary conditions. For this reason, we 
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revisit these representations in j66] . where we present a simplification, as well as a 
significant extension, of the results of [TJ |35] . 

1.2. Linear evolution PDEs with a second versus a third order x-derivative. 

Our definition of an effective characterization of A(k) and B{k), suggests that before 
attempting to solve a nonlinear boundary value problem, it is imperative that we have 
first constructed an effective solution of the underlying linearized problem. In this 
respect we note that linear evolution PDEs with second order derivatives are rather 
special, namely they can be solved by certain transforms which cannot be applied 
to linear evolution PDEs with third or higher order derivatives. In this sense, there 
exists a crucial difference between the linearized version of the NLS, namely of the 
equation 

iut + Uxx = 0, (1.1) 
and the linearized version of the mKdV, namely of the equation 



Ut + Uxxx = 0. 



(1.2) 



Let us consider the Dirichlet problem on the half-line for equations (1.1) and ( |1.2[ ), 
i.e. let us consider (1.1) and (1.2) in the domain 

< X < oo; < t < T, T positive constant, (1-3) 

with 

u{x,0) = uo{x), < X < oo; u{0,t) = go{t), 0<t<T, (1.4) 

where uo{x) and go{t) are given functions with appropriate smoothness and decay, 
satisfying uo{0) = go(0). It is well known that the classical sine transform can be used 



for solving equation (1.1) with the conditions (1.4). We recall that this transform can 



be obtained from the application of the Fourier transform to an odd function. Thus, 
the application of the sine transform is equivalent to solving the above problem via an 
odd extension of u from the half-line to the full line. Similarly, it is possible to solve 



equation (1.1) on the half-line with Neumann boundary conditions, either by using 



the classical cosine-transform, or equivalently via an even extension of u from the 
half- line to the full line. The solution of the Robin boundary value problem can also 
be constructed by either using an appropriate x-transform, or via a certain extension 
of u from the half-line to the full line. 



In contrast to equation ( 1.1 ), there does not exist a classical x-transform for solving 



equation (1.2) on the half-line. 



The unified method provides an elegant and effective solution for both equations 



(1.1) and (1.2). Indeed, define the following transforms of the given functions uo(x) 
and go{t), 



uo{k) 
Then, 



-ikx 



uo{x)dx, Im k < 0; go{k) 



e''^go{T)dT, k£C. 



u{x, t) 



1 

2^ 



^ikx—iuj(k)t 



UQ(k)dk 



1 

2^ 



^ikx—iu){k)t 



dD+ 



g{k)dk, 



where the (x,t)-domain is defined in (1.3), u}{k) is the dispersion relation, dD 
oriented boundary of the domain D 



:i.5) 



is the 



where D'^ is defined by 
D+ = {Im /c > 0, Im uj(k) > 0}, 

is to the left of the increasing orientation, see Figure 



with the assumption that D 
[l| and g{k) can be explicitly computed in terms of uq and (Jq: 
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Im k 



Im k 



Re k 
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kJl 







Re k 



Figure 1. The contour dD^ in the complex k-plane for equation (1.1), figure (a), 



and equation (1.2), figure (6). 



Linearized NLS 



uj{k) = k^ 



g{k) = uo{-k) - 2kgo{ik'^ 



Im A; > 0. 



(1.6) 



Linearized mKdV 



-ik-" 



a 



'Hit 

e 3 , Im K > 0. 



ijj{k) = —k^, g{k) = — ano(aA;)— Q^no(a^A;)+3A;^^o( 

By deforming the contour of Figure [T}x from the positive imaginary axis to the 
negative real axis, it is possible to recover from equation ( |1.5[ ) the classical solution 
obtained via the sine-transform (or equivalently via an odd extension): 



2 f°° 
u{x, t) = — ( 
Jo 



sin kx 



sm{k^)uo{^)d^ + ik I e^^ ^ gQ{T)dT 



dk. (1.7) 



However, as noted earlier, the right-hand side of equation (1.7) is not uniformly con- 
vergent at a; = 0. Indeed, if it were uniformly convergent then one could take the 
limit x — )• inside the integral and then u(0, t) would vanish, whereas we know that 
n(0,t) = go{t)- This lack of uniform convergence is a significant disadvantage of 
the solution obtained via a classical transform method for any boundary value prob- 
lem. This fact, in addition to rendering such representations unsuitable for numerical 
computations, also makes it difficult to prove existence without using PDF tech- 
niques. Indeed, taking into consideration that the construction of the solution via 
any transform method assumes existence of solutions, one must verify that the so- 
lution constructed under this a priori assumption satisfies the given PDF and the 
given initial and boundary conditions. An important advantage of the new method 
is that in addition to yielding an effective representation for evolution PDFs with 
spatial derivatives of arbitrary order, it also yields integral representations which are 
uniformly convergent on the boundary. This makes it straightforward to verify that 
the integral representation indeed solves the given initial-boundary value problem. 



For example, substituting (1.6) into (1.5) and evaluating the resulting equation at 
we find 



0, 



u{0,t) 



1 

2^ 



U{)[k)dk 



1 

2^ 



-ik^t 



dD+ 



Uo{—k)dk 



+ 



2ke 



-ik^t 



ik' 



''go{T)dT dk. 



(1.8) 



1 

By deforming the contour of integration in the second integral on the right-hand side 



T 



of (1.8) from the positive imaginary axis to the negative real axis and then replacing 
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k with —k in the resulting integral, it follows that the first two integrals on the 



right-hand side of (1.8) cancel; by making the change of variables k'^ = I in the third 



integral on the right-hand side of (1.8) and then using the integral formula of the 



classical Fourier transform, it follows that u{0,t) = go{t). 

1.3. Organization of the paper. In section [2] we review the main steps of the 
unified method; proofs of the basic results can be found in |3T1I35]. In sections |3] and 
|4]we derive effective characterizations of the spectral functions A{k) and B{k) for the 
Dirichlet and Neumann boundary value problems for the NLS and mKdV equations, 
respectively. In section [5] we utilize these characterizations to give a perturbative 
construction of the Dirichlet to Neumann map. In section |6] the above results are 
discussed further. 

1.4. Basic assumptions and notations. 

• We assume that q{x,t) vanishes sufficiently fast for all t as x ^ oo. Further- 
more, we assume that the given initial conditions satisfy 

qo{x) := q{x,0) G Li(M+) n L2(M+). (1.9) 

• 90i 91: 92 will denote q, Qx, and qxx evaluated at j; = 0, i.e. 

go{t) = q{0,t), gi{t) = qx{0,t), <?2(t) = fe(0, t), < t < T. (1.10) 

• Let Ai and A2 denote the two column vectors of the 2x2 matrix A. The 
notation ^(A;), k G {Di,D2), means that for Ai{k), k G Di, and for A2{k), 
k £ D2. 

2. Preliminaries 

Several important integrable nonlinear PDEs admit the following Lax pair |60j 
formulation: 

— {x,t,k) + ifi{k)a3fj,{x,t,k) = Q{x,t,k)fi{x,t,k), (2.1a) 
ox 

dfx ~ 

— {x,t,k) + if2{k)a3fi{x,t,k) = Q{x,t,k)fi{x,t,k), a;,tGM, k £ C, (2.1b) 

where /i is a 2 x 2 matrix-valued function, {fi{k), /2(fc)} are given analytic functions 
of k, and the 2x2 matrix- valued functions Q and Q are given analytic function of k, 
of q{x,t), of q{x,t), and of the x-derivatives of these functions. The action of CJ3 on 
a 2 X 2-matrix A is defined by 

<73^= [^73,^], a3 = diag(l,-l), (2.2a) 

and hence 

g*3x^ = e-'^x^g-^a^ = ( _Aii ^^"l^isA _ .2.2b) 

A21 A22 J 

Example 1 For the NLS equation, 
^| + 0-2AkPg = O, A = ±l, (2.3a) 

h{k) = k, f2{k) = 2k\ Qi^^t) = [xq{X:t) ^^^0*0' ^^'^^^ 
g(x, t) = 2kQ{x, t) + Q(i)(x, t), Q^^\x, t) = -i{Qx{x, t) + \\q{x, t)|^)a3. (2.3c) 
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Example 2 For the niKdV equation, 

q real, A = ±1 

h{k) = k, f^{k) = Ak\ Q{x,t)-- 



ot ox-^ ox 



q{x,t) 
\q{x,t) 



(2.4a) 
(2.4b) 



Q{x, t) = 2Q{x, tf - Q^^(a;, t) - 2ik{Q{x, tf + Q^{x, t))a3 + Ak'^Qix, t). (2.4c) 



2.1. The basic eigenfunctions. Using the identity (2.2b), it follows that equations 
(2.1) can be written in the form 

d 



dx 
d 

Ft' 



j{h{k)x+Mk)t)a3 , 



e-— ■ — -^/u(x, t, k) = eH^i W^+/2ifcjt)a3Q(2;, t, /t)^(x, t, k), 



AMk)x+f2(k)t)&3 



k £ 



Hence, 



d 



where the differential form W is defined by 

W{x, t,k) = ( Q{x, t, k)dx + Q{x, t, k)dt ) /x(a;, t, k). 



(2.5a) 
(2.5b) 



G C, 
(2.6a) 

(2.6b) 



Equation (2.6a) implies that if (x, t) G fi, where C is a simply connected 
domain, then the following equation characterizes a function /ij(x, t, k) which satisfies 
both ([2la|) and (I^Jbl): 



Hj{x,t, k) = I + 



(x,t) 



{Xj,tj) 



-i(/l(fc)(^-e)-3-i/2(fc)(t-r)a3p^^.(^^^^^)^ (x,t),(Xj,tj) G 17, 

(2.7) 



where is the differential form defined in (2.6b) with /i replaced by and {xj,tj) 
is a fixed point in fi. 

The implementation of the unified method relies on the fact that if $7 is a convex 
polygon, then by choosing {xj,tj) as the collection of the vertices of this polygon, it 
follows that the collection of fij defines a sectionally analytic function in the complex 
fc-plane 

In the case of the half-line, is given by 



= {0 < X < oo, < t < T}, 
thus the relevant polygon is convex and the three corners are the points 

(0,T), (0,0), (oo,t). 



(2.^ 



We refer to these corners as {1, 2, 3}, respectively. Let {fii, ^2, fJ's} be the eigenfunc- 
tions corresponding to these corners. These eigenfunctions are independent of the 
contours from {xj,tj) to (x,t). Typical contours are shown in figure [sj 
For fJ-i,H2, /^3, the following inequalities are valid respectively: 



x-^ > 
t-T <0 



X - ^ > 
t-r > 



x-^ < 
t-r = 
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w 



• 3 



X 



Figure 2. The domain Q with the corners at (0, T), (0,0), and (oo,t). 




T 



{x,t) 



(1) (2) 
Figure 3. The contours used to define {/ijlf- 



(3) 



(a) 



Da 



-Re k 




Figure 4. The domains {Dj}f for (a) NLS and (b) mKdV. 



Re k 



Hence, the expressions that appear in the first and second columns of (2.7) are analytic 
and bounded in the following domains of the complex /c-plane: 



(Im /i > n Im /2 < 0, Im /i < n Im /s > 0) =: {D2, D3), 
(Im /i > n Im /2 > 0, Im /i < n Im /s < 0) =: iDi,D^), 
(Im/i <0, Im/i >0) =: (C-,C+). 



fJ'2 

For NLS and mKdV the domains {Dj}f are shown in figure |4| For these equations, 
fi{k) = k so that C~ and C"^ denote the lower and upper halves of the complex 
/c-plane, respectively. 

2.2. The spectral functions. The crucial advantage of requiring that the functions 
solve equation (2.7) is that any two such functions are related via a matrix which 
has explicit exponential [x, t)-dependence of the form exp[—ifi{k)xa3—if2{k)ta3]p(k), 
where the function p{k) can be computed by evaluating the relevant relation at any 
convenient point in the domain $7. In particular, using that //2(0,0,/c) = /, we find 
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the following relations: 

/X3(x,t,fc) = ^2(x,t,fc)e-^(^iW^-+^2('=W^«s(A:), k E (C-,C+), 
t, k) = fi2{x, t, A;)e-^(^^ k G {D^, D^), 

where the spectral functions s{k) and S{k) are defined by 

s{k) = ^{0,k), A:E(C-,C+), 
S^^{k) = e'f^^^^^{T, k), keC, 



(2.9a) 
(2.9b) 



(2.10) 
(2.11) 



and <I> denote the functions obtained by evaluating /xa and /i2 at t = and x = 
respectively, i.e. 

^{x,k) = iJ3{x,0,k), A:e(C",C+); ^{t,k) = fi2{0,t,k), keC 

Hence, ^ and satisfy the x and t-parts of the associated Lax pair evaluated at t = 
and X = respectively: 



and 



+ i/i<T3^' = Qo^', Qo{x,k) = Q{x,{),k), 
lim^' = /, 0<x<oo, fcG(C~,C+), 



^t + if2^z^ = Qo^, Qoit,k) = Q{0,t,k), 
^{0,t) = I, 0<t<T, keC. 



(2.12) 



(2.13) 



Equations (2.12) and (2.13) are equivalent to the following linear Volterra integral 
equations: 



^{x,k) = I 



{Qo^)i(.,k)d^, 0<x<oo, k£{C~,i 



(2.14) 



and 



$(t, k)=I + [ e-'^2(fc)(t-r)a3(g^^)(^^ < t < T, k^C. (2.15) 

Jo 

2.3. The global relation. Let F(t, k) denote the eigenfunction obtained from the 
evaluation of /X3 at x = 0, i.e. 

F{t,k) = ti3{0,t,k), A;G(C-,C+). 



The eigenfunctions and F satisfy the same differential equation (2.1b), thus they 
are simply related: 

^{t,k)e-'^^^'''^^^^s{k) = F{t,k), 0<t<T, ke{C-,C+). (2.16) 
Evaluating this equation at t = T we obtain the global relation 

S-^{k)s{k) = e'^^^^'F{T, k), ke (C^, C+). (2.17) 
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2.4. Notations. The superscripts (+, — , 1,2,3,4) will denote the domain of analyt- 
icity of the relevant vector eigenfunction: 

^l = (^f^M?^)' ^2 = (4^\m2^^) , H3 = {fJ-s , fit) ■ (2-18) 

For NLS and mKdV, the matrices Q and Q possess certain symmetries, which in turn 
imply certain symmetries for the eigenfunctions and for the spectral functions: 

^t, k) = ( ^Sl , Sik) = ( ^(^)'l , keC, (2.19) 



^A^i(x,/c) ^2{x,k)j' \\b{k) a{k)^ 

We will henceforth assume that these symmetries are present. 



(2.20) 



2.5. The RH problem. Equations (2.9) together with the estimate 

/ij. = / + oQ^, A:^oo, i = 1,2,3, (2.21) 

for k in the domain in which /ij is bounded and analytic, can be used to formulate a 
RH problem. Indeed, equations ( |2.9| ) can be written in the form 

M^{x,t,k) = M+ix,t,k)J{x,t,k), k€D+r\D^, (2.22) 

where 

D+ = DiU Ds, D- =D2U Di, 
and the matrices M_, M^, J are defined as follows: 

M,= (|i,4), keD,, M_=(g,4), k,D2; 

M+ = f/X3 , d^] , keDs; M_ = ( 1^3,4^] , kGD^; 
V d{k)J V 

'Ji, kGDinD2, 
J2 = hJi^Ji. keD2nD3, 

J3, kGDsD D4, 

J4, kGDinDi; 



J(x, t, /c) 



with 



a{/c) 
1 

a{k)a(k) 



d{k) = a{k)A{k) - Xb{k)B{k); e{x, t, k) = fi{k)x + f2{k)t. 

2.6. The main rigorous result. The above derivation is based on the assumption 
that there exists a solution of the given nonlinear PDE with appropriate smoothness 
and decay. The rigorous justification of the unified transform method is based on the 
following important result, first proved in [35] in connection with the NLS equation: 
(i) Given initial condition qo{x) G Li(M"'") H L2(M"'"), define the spectral functions 
{a{k),b{k)} hy 

a{k) = ^-2(0, k), b{k) = ^-1(0, k), ke C+, 
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where "^{xjk) is the solution of ( 2.12[ ), or equivalently of (2.14) and QQ{x,k) is 
uniquely defined in terms of qo{x) (with q{x,0) replaced by qoix))- (ii) Given a 
subset of the boundary values 

go{t) = q{0,t), gi{t) = q^{0,t), 52(i) = fe(0, t), 

assume that it is possible to construct the spectral functions {A{k), B{k)}, 

A{k) = ^2{t,k), B{k) = -e^'f^^^i{T, k), keC, (2.23) 



un der th e following requirements: (a) ^{t,k) is a solution of (2.13), or equivalently 
of (2.15), where Qo{t,k) is defined in terms of gj, j = 0,1,2. (6) $ satisfies the 
restriction 



^t,k)e-'f^^''^^^^s{k) = F{t,k), 0<t<T, A;G(C-,C+), 



(2.24) 



where F{t,k) is some function whose first and second columns are analytic in k in 
C~ and C"*" and 



F{t,k) = 1 + 



CXD, 



k G (C-,C+). 



(iii) Define the RH problem with the jump condition (2.22) in terms of the spectral 
functions {a{k) , b{k) , A{k) , B{k)}. Then, the unique solution of this RH problem can 
be used to compute a function q(x, t) which solves the given nonlinear PDE and 
furthermore. 



q{x,0) = qoix), 



< X < oo; 



d^q 



{0,t)=gjit), 



0,1,2, 



< t < T. 

(2.25) 



The above theorem shows that the rigorous justification of the unified method 
reduces to the following problem: Given {a{k),b{k)} and a subset of {go, gi, g2}, 
characterize {A{k), B{k)} from the requirement that ^{t, k) satisfies (2.13) and ( 2.24[ ). 

It turns out that: (a) For linearizable boundary conditions, it is possible to express 
explicitly {A(k),B{k)} in terms of {a{k),b{k)} and the given boundary conditions. 
(b) For non-linearizable boundary conditions, it is possible to express {A(k),B{k)} 
in terms of a system of two nonlinear equations which can be iteratively solved for 
'small' boundary conditions. 

Remark 2.1. For the NLS, there exists one unknown boundary value. Thus, one 



might expect that since the global relation (2.24) provides one equation connecting 



this unknown function with the given initial and boundary conditions, in this case 
it is possible, by utilizing the global relation, to characterize the unknown boundary 
value. However, for the mKdV, there exist two unknown boundary values. In this 
case, the solution of the associated linear problem suggests that it is impossible to 
solve this problem, unless one uses the transformation which leaves f2{k) invariant in 
order to obtain an additional equation from the global relation. 

Remark 2.2. It is important to note that the rigorous result mentioned earlier does 
not require the knowledge of the explicit form of F; it only requires the existence of a 
function F with specific analyticity properties. This suggests that the most efficient 
approach of characterizing {A{k),B(k)} is to actually eliminate F. This is precisely 
the philosophy used in the linear limit for the derivation of the generalized Dirichlet 
to Neumann map and it is also the philosophy used in ^ [351 US] ■ 
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3. The NLS equation 

We consider the Dirichlet and Neumann boundary value problems for the NLS 
equation (2.3a) posed on the half-line. 

3.1. The global relation. Define c{t,k) by 

c{t,k) = ^{F{t,k))i2, (3.1) 

where F{t, k) is some function whose first and second columns are analytic and 
bounded in and respectively and F{t,k) = / + 0{l/k) as /c — )• cxd with 
k S (C~,C^). We can write the (12) entry of the global relation (2.17) as 



b{k)- 



c{t, k) = $i(t, k) + ^^2(t, k)e 
a[k) 



'2if2{k)t 



Im/c > 0. 



(3.2) 



The function c(t, k) is analytic and bounded in Im k > away from the possible zeros 
of a{k). The functions a{k) and b{k) are analytic and bounded in D1UD2, while A(k) 
and B{k) are entire functions which are bounded in Di U D3. The functions ^i{t, k) 
and ^2it,k) are entire functions which are bounded for k S D2 U D4. Furthermore, 
the second column vector of equation (2.15) yields 
ft 



'Aik^{t~T) 







(Qo)ii^i + (Qo)i2$2 {r,k)dT, 



^2{t,k) = 1 + 
Hence, 



t r- 



(Qo)2i^i + (<9o)22^2 (r, k)dT, 



0<t<T, ke 



(3.3) 



ii(i,A:) 
where 



((3o)ii«'i + ((5o)i2^>2 {T,k)dT, 



((3o)2iii + ((§0)22^2 (r,A;)dr, 



^j{t,k) = <!>j{t,k)e 



iik'^t 



j = 1, 2. 



(3.4a) 

< t < T, A; G C, 
(3.4b) 

(3.5) 



Therefore, $j(t, k), j = 1, 2, is bounded for k G DiVJD^ and ^j{t, k)e '^^'^ * is bounded 
for k G D2UD4. Equation (3.2) shows that whereas <I>i and ^26"^*^^ * are bounded in 
D2 U D4, the combination appearing on the right-hand side of (3.2) is also bounded 
in Di. 



3.2. Asymptotics. Integration by parts in (3.3) shows that 



^i{t,k) 



+ 



k^ 



+ 



+ 



-Aik'^t 



k 



A; — )• 00, k £ D2U D4, 



<^2it,k) = l+ 2 
where 



k ^ 



9o{t) 
2i ' 



(0,0) 



A(e,T), 



k ^ 00, k e D2LI Z?4, 

giit) igoit) /■(°'*) 



(3.6) 



/ A(e,r) 

^(0,0) 
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and the closed one-form A is defined by 

^{x,t) = ^ [-i\q{x,t)\'^dx + {q{x,t)qx{x,t) - q{x,t)q.j:{x,t))dt] . 
In particular, we find the following expressions for the boundary values: 

go{t) = 2i^']'\t), (3.7a) 

gi (t) = ) (t) + 2^50^^'^ (t), 0<t<T. (3.7b) 
We will also need the asymptotics of c. 



Lemma 3.1. The global relation (3.2) implies that the large k behavior of c{t,k) 
satisfies 

c{t,k) = ^^^ + ^\P- + o(^), k^oo, keDi. (3.8) 



k k'^ \k^. 

Proof. See appendix [Bj □ 

3.3. The Dirichlet and Neumann problems. The following theorem expresses 
the spectral functions A{k) and B{k) in terms of the prescribed initial and boundary 
data via the solution of a system of nonlinear integral equations. This result was 
already derived in [33] using the GLM representations for $i and ^>2- However, the 
derivation here does not require the use of the GLM representations; the result is 
instead deduced directly from the asymptotics of the eigenfunctions as A; — )• oo. 

We let dDj, j = 1, ... ,4, denote the boundary of the j'th quadrant Dj, oriented 
so that Dj lies to the left of dDj. 

Theorem 3.2. Let T < oo. Let qo{x), x > 0, be a function of Schwartz class. For 
the Dirichlet problem it is assumed that the function go{t) , < t < T , has sufficient 
smoothness and is compatible with qo{x) at x = t = 0. Similarly, for the Neumann 
problem it is assumed that the function gi{t), < t < T, has sufficient smoothness 
and is compatible with qo{x) at x = t = 0. Suppose that a{k) has a finite (possibly 
empty) set of simple zeros, which are denoted by {kj^i ; assume that no zeros occur 
on the boundaries of Di and D2 . 



Then the spectral functions A{k) and B{k) associated with the NLS equation (2.3a) 
are given by 



A{k) = $2(T, k) B(k) = -«>i(T, k)e^''' ^, (3.9) 

where the complex-valued functions ^i{t,k) and ^2{t,k) satisfy the following system 
of nonlinear integral equations: 
ft 



^i{t,k)= [ e'"'"'-''-'^[-iX\go\^^i + i2kgo + igi)^2]it',k)dt', (3.10) 
Jo 

<^2{t, k) = l + \ [ [{2kgo - %)$i + i|ffoP$2] {t', k)dt', 0<t<T, keC. 
Jo 



(a) For the Dirichlet problem, the unknown Neumann boundary value gi{t) is 
given by 

ai{t) = - [ {kxi{t, k) + igo{t))dk + f X2{t, k)dk 

JdDz ^ JdD-i 



JdD^ a{-k) 
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kjeDi 



a{kj) 



^2{t,kj), 0<t<T, 



(3.11) 



where {Xj}i denote the odd combinations formed from i.e. 
Xj{t,k) = ^j{t,k) - ^j{t,-k), i = l,2, 0<t<r, keC. (3.12) 
(6) For the Neumann problem, the unknown Dirichlet boundary value go{t) is 



given by 



1 



TT 



go{t) = - / xi{t,k)dk + - 



dD3 



TT 



e-^"'''^^<S>2it,-k)dk (3.13) 



dD3 



a{—k) 



ik'it 



e 3 



d{kj 



^2{t,kj 



where {Xj}i denote the even combinations formed from i.e. 
Xj{t,k) = ^j{t,k) + <^j{t,-k), j = l,2, 0<t<T, keC. (3.14) 



Proof. Equations (3.9) and (3.10) follow from the definition of $i and <I>2 and from 
(I2l5l). 



(a) In order to derive (3.11) we note that equation (3.7b) expresses gi in terms of 



#2^^ and Furthermore, equations (3.6) and Cauchy's theorem imply 



dD2 



[<l>2{t,k) - l]dk 



[^2{t,k) - l]dk 



(3.15) 



and 



-f*f'(*) 



dD2 



k^i{t,k) 



9o{t) 
2i 



dk 



k^i{t,k) 



9o{t) 
2i 



dk, (3.16) 



where we have used equation (3.7a) to write ^^^^ in terms of go{t). 

If, instead of equation (3.11), we use equation (3.7b) with <I>2^^ and ^^^^ defined 
by (3.15) and (3.16), we do not obtain an effective characterization. The latter 
representation requires the appearance of xi X2- In this respect we note that 



Z7T'^^i\t) 



Similarly, 
iTr<^f\t) 



+ 

dD2 JdD4, 



[<^>2it,k) - l]dk 



+ 

dDz JdDi 



[$2(t,fc) - l]dk 



[ [^2{t,k)-l]dk- f [^2{t,-k)-l]dk= f X2it,k)dk. (3.17) 
JdDs JdDs JdDs 



+ 



'dD3 

k<^>i{t,k) 
k<^>i{t,k) 



gojt) 

2i 

9o{t) 
2i 



dk 

dk + 2 



dDi 



k<^>i{t,k) 



gojt) 

2i 



dk 



[kxiit,k) + igo{t)]dk + 2 



dDi 



k<^i{t,k) 



9o{t) 



2i 



dk. 



(3.18) 



The last step involves using the global relation (3.2) in order to compute the second 

go{ty 



term on the right-hand side of (3.18): 

goity 



dDi 



k^i{t,k) 



2i 



dk = 2 



dDi 



kc{t, k) 



2i 



dk 
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2 / ^ci,2(t,fc)e-2^/2W*dA;. 



(3.19) 



Using the asymptotics (3.8) of c{t, k) and Cauchy's theorem to compute the first term 
on the right-hand side of (3.19) and using the transformation k — >• —k in the second 

kb{-k) 



term on the right-hand side of (3.19), we find 

go{t) 



k^i{t,k) 



2i 



dk 



a{-k) 



^2it,-k)e-^'^'^ 



a{kj) 



^2{t,kj)e-^'^'K 



(3.20) 



(6) In order to derive (3.13 
Furthermore, equation ( 



Equations (3.17), (3.18), and (3.20) together with (3.7b) yield (3.11). 



we note that equation (|3.7a|) expresses qq in terms of 

(3.21) 



3.6) and Cauchy's theorem imply 
^i{t,k)dk 



aD2 



dD4 



^i{t,k)dk. 



Thus, 



+ / ]<^i{t,k)dk 
dDz JdDi, 



<^i{t,k)dk + 2 / <^i{t,k)dk 

dDi/ JdDi 



Xi{t,k)dk + 2 



dDi 



^i{t,k)dk. 



(3.22) 



The last step involves using the global relation (3.2) and to compute the second term 



on the right-hand side of (3.22): 

2 / ^i{t,k)dk = 2 c{t,k)dk 
JdDi JdDi 



dD 



m 

^a{k) 



^2{t,k)e-'^'f^^^^^dk. (3.23) 



Using the asymptotics (3.8) of c{t, k) and Cauchy's theorem to compute the first term 
on the right-hand side of (3.23) and employing the transformation A; — )• —km. the 
second term on the right-hand side of ( 3.23| ), we find 

h{-k) 



2 I <^i{t,k)dk 
IdDi 



iTT^'l\t) + 2 



dD3 



a{-k) 



^2{t,-k)e- 



fc,e-Di 



MM 

a{kj) 



^2{t~kj)e-^'^'' . 



(3.24) 



Equations (3.22) and (3.24) together with (3.7a) yield (3.13) 



□ 



Remark 3.3. The Dirichlet and Neumann problems for the NLS equation on the 
half-line can now be solved as follows: If a{k) has no zeros, the functions A{k) and 



B{k) enter the formulation of the RH problem of subsection 2.5 only in the combina- 



tion A{k) / B{k) for k E dD^. For the Dirichlet problem, substituting the expression 



(3.11) for gi{t) into (3.10) yields a system of quadratically nonlinear integral equa- 



tions involving the functions $i(t,fc) and ^>2(ii^)) 0<t<r, /cSMU zM. For the 



Neumann problem, substituting the expression (3.13) for gQ{t) into (3.10) yields a sys- 



tem of quadratically nonlinear integral equations for <I>i(t, k) and $2(^5 k). Assuming 
that these systems have unique solutions, A{k)/ B{k), k € dD-^, can be determined 
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from (3.9). In fact, we will show in the following subsection that these nonlinear 
systems provide effective characterizations of A{k) and B{k) in the sense defined in 
the introduction. In particular, the systems can be solved recursively to all orders in 
a perturbative scheme. If a{k) has zeros, the sum of residues also needs to be taken 
into account. 



3.4. Effective characterization. Substituting into the system (3.10) the expan- 
sions 

^2 



2 



90 = £501 + £"502 + 
ho = ehoi + e^/io2 H , 



e^O, j = 1,2, 



gi = egu + e 512 H , 

hi = ehii + e^hi2 H 



where e > is a small parameter, we find that the terms of 0(1) give <I>io 
$20 ^ Ij the terms of 0(e) give $21 = and 

ft 



0(e) : ^ii{t,k) 
the terms of O(e^) give 

0(e2) : ^i2{t,k) = 



{2kgoi{t')+igu{t'))dt'; 



4ifc2(j'-t) 



{2kgo2it') + igi2{t'))dt'; 



(3.25) 
(3.26) 
(3.27) 
= and 

(3.28) 
(3.29a) 



0(e2): ^22{t,k) = x[ [{2kgoi-ign)^ii + i\goi\^]it',k)dt'. 
Jo 

On the other hand, expanding ( |3.11[ ) and assuming for simplicity that a{k) has no 
zeros, we find 



511 



2 

Tri 



{kxii{t,k) +igoi{t))dk 



4 
vri 



dD3 



-k)dk, 



(3.30) 



where xi = ^Xu + O(e^) and b = ebi + O(e^). 

The Dirichlet problem can now be solved perturbatively as follows: Let Xjj 3 = 1) 2, 
denote the odd combination in k formed from as defined in (3.12). The odd part 
of ( [31281 ) yields 

Xii{t,k) = Ak re4*'='(*'-*)5oi(t')rfi'- 



(3.31) 
Jo 

Given goi, we can use this equation to determine xii(^5^)- We can then compute 
gii from (3.30) and then $11 follows from (3.28). In the same way, we can use the 
odd part of (3.29a) to determine X12; we use xi2 to compute 512, and then ^>i2 and 
$22 follow from ( 3.29| ). This recursive scheme can be continued indefinitely. Indeed, 
suppose and gij have been determined for all < j < n — 1, for some n > 0. 

The terms in (3.10) of 0{e^) give 

ft 



0{e'') : ^in{t,k)= [ 
Jo 



e'^ik^{t'-t)(^2kgon{t') + igin{t'))dt' + lower order terms, 

(3.32a) 
(3.32b) 



0{e^) : ^2n = lower order terms, 

where 'lower order terms' denotes an expression involving known terms of lower order. 
Similarly, the terms of 0(e") of the integral representation (3.11) for gi give 



9ln 



7r« JdD. 



{kxin+igon)dk- 



4 

TTl 



ke '^'^ %n{-k)dk+\awex order terms. (3.33) 



dD-j, 
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The odd part of ( |3.32a ) yields 



Xlnit,k) 



Ak I e 

'0 



Aik'^(t'-t) 



gon{t')dt' + lower order terms. 



Substituting the solution xin of this equation into (3.33), we find gin and then <I>i„ 



and ^2n are found from (3.32). This shows that for the Dirichlet problem $i and 



<I>2 can be determined to all orders in a perturbative scheme by solving the nonlinear 
system of theorem |3.2| recursively. 



Similarly, for the Neumann problem, substituting the expression (3.13) for go{t) 
into (3.10) yields a system of quadratically nonlinear integral equations for the func- 



tions ^i{t,k) and ^2{t,k). This nonlinear system can be solved recursively to all 
orders perturbatively and therefore it provides an effective characterization of A{k) 
and B{k) for the Neumann problem. Indeed, letting 



Xi = eXii + eXi2 + 



the even part of (3.28) yields 

Xii{t,k) 



2i I e 
'o 



gu{t')dt', 



while (3.13) yields 
ffoi 



- / xii{t,k)dk + - I e 
^ JdD3 ^ JdDz 



bi{-k)dk. 



(3.34) 



(3.35) 



Since gi is known, (3.34) can be solved for xii) and then g^i can be determined from 



(3.35). After g^i has been found, $ii can be computed from (3.28). Extending this 



procedure to higher orders we find, just like in the case of the Dirichlet problem, a 
recursive scheme which can be used to determine $i and $2 to all orders. 



Remark 3.4. The function g^, in addition to (3.13), also admits the following alter- 
native representation: 

9o{t) = - [ Xiit,k)dk. (3.36) 

However, this alternative representation is not suitable for the effective solution of 
the Neumann problem. Indeed, for the Neumann problem, we can find xin from the 



terms of 0{e'^), so that the representation (3.13) can be used to find gon- However, 
the function xin remains unknown, so that the representation (3.36) cannot be used 
to find gon- 

4. The mKdV equation 



We consider the mKdV equation (2.4a) posed on the half-line 



4.1. Asymptotics. In this case, the $,'s admit the asymptotics 



$i(t,A;) 



(3), 



k 



k^ 



+ 



fc3 



+ 



+ 



Mt,k) = 1 
where 



cl>f (t) / 1 

k k^ U3 



(1) 



2i 



A, 



(2) 



(0,0) 



(0,t) 



(0,0) 



(1) 



00, A;G£)2U£>4, 



go 

2i' 



(4.1a) 
(4.1b) 



18 



(2) _ 91 Affo 



A. S. FOKAS AND J. LENELLS 
(0.*) ,n^ 1 



A, 



(0,0) 



(3) _ M „ 3,(2) , i ^ 3,(1) , A^o 
1 -^Uo<&2 +251^2 +^ 



r52 



and the closed one-forms A and A are defined by 



A = q^dx + {q^ — 2qqxx + 3Ag^)di, 



A 




In particular, we find the following expressions for the boundary values: 

(1) 



50 = 

.(3) 



92 = \9^ - ^i^r + 450^f + 2igi<^ 



.(1) 



We will also need the asymptotics of the function c(t, k) defined in (3.1 ). 



(4.2a) 
(4.2b) 
(4.2c) 



Lemma 4.1. The global relation (3.2) implies that the large k behavior of c{t,k) 
satisfies 

.,,„.£i^,iSM,£S!^,0(J,). ..oo,..«.. ,4.3) 
Proof. See appendix [B] □ 



4.2. The Dirichlet and Neumann problems. We can now derive effective char- 
acterizations of A{k) and B(k) for the Dirichlet (^o prescribed), the first Neumann 
{gi prescribed), and the second Neumann {g2 prescribed) problems. 

Define a by a = e^'^^^^ and let {XjjXj}Xj}i denote the following combinations 
formed from {<&j}i: 

Xj{t, k) = ^j{t, k) + a^j{t, ak) + a^(^j{t, a^k), j = 1, 2, 

Xj{t, k) = ^j{t, k) + a^(tj{t, ak) + a$j(f, a^k), j = 1, 2, 

Xj{t, k) = ^j{t, k) + $j(i, ak) + $j(t, a^k), j = 1, 2. 

Define i?(t, k) 



R{t,k) 



b{k)^2{t,k) 
a{k) ' 



< t < T, k£ 



Let Di = D[UD'{ where D[ = Din{Rek> 0} and D'{ = Din{Re A; < 0}. Similarly, 
let L>4 = U D'l with D'^ = D^n {Re A; > 0} and D'l = L>4 n {Re A; < 0}. 

Theorem 4.2. Let T < oo. Let qo{x), x > 0, be a function of Schwartz class. For 
the Dirichlet problem it is assumed that the function gQ{t) , < t < T , has sufficient 
smoothness and is compatible with qo{x) at x = t = 0. Similarly, for the first and 
second Neumann problems it is assumed that the functions gi{t) and g2{t) , < t < T , 
have sufficient smoothness and are compatible with qo{x) at x = t = 0, respectively. 
Suppose that a{k) has a finite (possibly empty) set of simple zeros, which are denoted 
by {kj}i ; assume that no zeros occur on the boundaries of Di and D2. 

Then the spectral functions A{k) and B{k) associated with the mKdV equation 
(2.4a\ ) are given by 



A{k) = $2(r, k) B{k) = -^>i(r, k)e 



(4.4) 



The f-dependence of R{t, k) will be suppressed below. 
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where the complex-valued functions ^i{t,k) and ^2{t,k) satisfy the following sysl 
of nonlinear integral equations: 
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,em 



,{t, k) = f eS^'^'^*'-*) [-2ik\gl<^^ + {2\gl + 4^:^50 + 2i%i - 52) $2] k)dt' , 
Jo 



(4.5a) 



$2(i, k) = l + \ [ [{2Xgl + 4pgo - 2ikgi - 92) $1 + 2ikgl(t2] {t\ k)dt' , 
Jo 

0<t<T, k€C. (4.5b) 

(a) For the Dirichlet problem, the unknown Neumann boundary values gi{t) and 
92{t) are given by 



g,^t)=^ [ X2it,k)dk+- 



TT JdDa 

-—I ke-^'^'^ \{a^ - a)R{ak) + {a - a^)R{a^k)\ dk 

JdDz 



kxi{t, k) 



3go{t) 



2i 



dk 



and 



+ 4^1 -a') Yl +(!-") E W 

3kgo{t) 



-8ikH 



Res R{k) 



(4.6a) 



92{t) = Xglit) 



+ 



- / k'^e-^'''^*Ul-a)R{ak) + {l-a'^)R{a'^k)]dk 



k\i{t,k) 



2i 



dk 



-8^(1 -a) J2 +(!-«') E V. 



251 (i) 



^e-^^'^'ResRik) 



+ ^/ kUt,k)dk + 



X2{t,k)dk. 



(4.6b) 



(6) For the first Neumann problem, the unknown boundary values go{t) and g2{t) 
are given by 

5o(i) = -/ xi{t,k)dk-- I e-^^^^^{a - a^)R{ak) + [a^ - a)R{a^k)\dk 

+ 2i|(l-a) +{l-a^) Y |e"*'''^'*Rcsi?(A;) (4.7a) 



kjeD'i 



and 



g2{t) = Xglit) - - / (k^xi{t,k) - - / lxi{t,l)dl) dk 

[ k'^e-^'^''* [(1 - a'^)R{ak) + (1 - a)i?(a2it)] dA; 



71" JaD3 
8z 



2-8ik-H 



Res R{k) 



+ 



.{(l-a^) 5^ +(l-a) Y ]k: 

[ kx2{t,k)dk+'^^ [ X2{t,k)dk. (4.7b) 

JdDa Jqd 



(*) 
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(c) For the second Neumann problem, the unknown boundary values go{t) and 
gi{t) are given by 



go 



{t) = - [ xi{t,k)dk-^ [ e-^''''^{a-l)R{ak) + {a^ -l)R{a^k)]dk 
^ JdD3 JdDa 

+ 2i|(l-a2) +(l-a) Yl ]e~^'''^'ResR{k) (4i 



and 



250 (t) 



vr 
2 

vri 



X2{t, k)dk + 



2 



dD-j, 



kxi{t, k)dk 



f ke~^''''' [{a^ - l)R{ak) + (q - l)R{a^k)] dk 



+ 4|(l-a) +(!-«') E y. 



e-'^^'^'ResRik). 



(4.8b) 



Proof, (a) In order to derive (4.6a) we not e tha t equation (4.2b) expresses gi in terms 
of <&2^^ and Furthermore, equations (4.1) and Cauchy's theorem imply 



2i7r (1) 



-—^''{t) = 2 / [<^2{t,k)-l]dk 



dD2 



and 



Thus, 



2m ^ (2) 



^'{>{t) = 2 



dD2 



k^i{t,k) 



gojt) 

2i 



dk 



dD4 



[<^2{t,k) - l]dk 



k^i{t,k) 



gojt) 

2i 



dk. 



ivrci>«(t) 



+ 

dD2 JdDiJ 



][^2{t,k) - l]dk 



+ 

dDz JdDi 



[^2{t,k) - l]dk 



[ [^2{t,k)-l]dk + a [ [^2{t,ak)-l]dk + a^ j [^2{t,a^k) - l]dk 

JdDi JdDs JdD-s 



X2{t,k)dk. 



(4.9) 



Similarly, 



ivr$f)(t) = ( / + / 

KJdDs JdDi 



k^i{t,k) 



gojt) 

2i 



dk 



kxi{t,k) 



dD[ JdD'{ 

3go(0 
2i 



k^i{t,k) 



dk + I{t), 



go{t) 



2i 



dk + I{t) 



(4.10) 



where I{t) is defined by 



+(l-a) 



dD' 



k^l{t,k) 



gojt) 

2i 



dk. 



The last step involves using the global relation (3.2) to compute I{t): 



I{t) 



I -a' 



+{l-a] 



dD'{ 



kc{t, k) 



gojt) 

2i 



dk 
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(l-a^) / +(l-a) 



dD'l 



(4.11) 



Using the asymptotics (4.3 ) of c {t, k) and Cauchy's theorem to compute the first term 
on the right-hand side of (4.11) and using the transformations k ^ ak and k — >• a'^k 



in the second term on the right-hand side of (4.11), we find 

I{t) = - ivr^P^ (t) - / ke-^''''^ [(a^ _ a)R{ak) + (a - a^)R{a^k)\ dk 

JdD-i 



27ri<j(l-a2) +(l-a) ^ \kje-^'''j'Res R{k). 



(4.12) 



Equations (4.10) and (|4.12|) imply 
1 



c|>f)(t) 



2TTi 



dD3 



kxi{t,k) 



2i 



dk 



1 



ke-^'^''^[{a^ - a)R{ak) + (q - a^)R{a^k)\dk 



dD3 



Res R{k). 



This equation togethe r wit h (4.2b) and (4.9 ) yie lds (4.6a). 



(3) (2) 

In order to derive (4.6b), we note that (4.2c) expresses §2 in terms oi ^\ , ^2 ' 



and $2^^ Equation (4.6b) follows from the expression (4.9) for <1>2^^ and the following 
formulas: 



(2) 
2 

(3) 



1 

1 

27ri 



kx2dk, 



(4.13a) 



k'^Xi 1 dk 



— [ k^e-^''''*\{l-a)R{ak) + {l-a'^)R{a^k)]dk 



2i 



+ /(!_«) ^ +(l-a2) ^ 



^e-^^'^'ResRik). 

kj 



(4.13b) 



(6) In order to derive the representations (4.7) relevant for the first Neumann 
problem, we use (4.2) together with (4.9), ( 4.13a| ), and the following formulas: 

^ ^ Xidk-— [ e-^''''*\(a-a^)R(ak) + (a'^ -a)R(a'^k)]dk 



(1) 



27ri 



dD-i 



|(l-a) +(!-«') Yl |e"^"'^'*Resi?(A:) 



(2) 



k,eD[ 
kxidk, 



$1' 



(3) 



2-Ki 



dDi 



(4.14a) 
(4.14b) 



k^Xi - ^) dk 



— [ k'^e~^'''^^\(l-a^)R(ak) + (l-a)R(a'^k)]dk 
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Re k 



Figure 5. The contour T is defined in such a way that aT and a^T pass above the 
zeros of a{k) (these zeros are indicated by x's in the figure). 



Ul-a^) +(l-a) Yl ]kje-^''''PResRik). (4.14c) 



(c) In order to derive the representations (4.8) relevant for the second Neumann 
problem, we use (4.2) together with (4.9) and the following formulas: 

1 



(1) 



2-Ki 



Xidk 



1 

2TTi 



-Sik'H 



dD3 



[{a - l)R{ak) + (a^ - l)R{a^k)\dk 



$1- 



(2) 



+ Ul-a^) Y +(!-«) |e"^*''^'*Resi?(fc), (4.15a) 
— f ke-^'''''[{a^ -l)R{ak) + {a-l)R{a^k)]dk 



1 

27H 



kxidk 



'dDz 
2\ 



+ {{l-a) Y +(!-«') Y 

kjGDi k,eD'(^ 



-Sikh 



Res R{k). 



(4.15b) 



The proofs of the formulas (4.13 )-(4.15 ) rely on arguments similar to those used in 
the proof of (4.6a). □ 



Re mar k 4.3 . The sums over the zeros of a{k) that appear on the right-hand sides 
of (4.6)-(4.8) can be absorbed into the integrals involving R{ak) and R{a^k) by 



replacing the integration contour dD^, for these integrals by F, where F is a contour 
obtained by deforming dD^, in such a way that aT and a^T pass above all the zeros 
of a{k) in D'^ and D", respectively, see figure 



4.3. Effective characterization. The nonlinear system for and ^>2 obtained 



by substituting the expressions (4.6) for gi and g2 into (4.5) provides an effective 



characterization for A{k) and B[k) for the Dirichlet problem for mKdV. Similarly, 
substituting the representations (4.7) and (4.8) into (4.5) yields an effective charac- 



terization of A{k) and B{k) for the first and second Neumann problems respectively. 
Indeed, substituting into the system (4.5) the expansions 



+ e^ii + e^j2 + • • • 



j = l,2, 
0,1,2, 
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where e > is a small parameter, we find that the terms of 0(1) give $10 = and 
$20 = 1) while the terms of 0(e) give ^21 = and 

ft 



0(e) : '^u{t,k) 



„8ifc3(t'-t) 







{4k'^goi + 2ikgn - 521) dt'. 



The symmetric combinations of equation ( |4.16[ ) are 



Xii{t,k) = l2k^ 



Xii{t,k) = 6ik / e^"^''-''-'^gn{t')dt', 



Xii{t,k) 



(4.16) 

(4.17a) 
(4.17b) 
(4.17c) 



The Dirichlet problem can now be solved perturbatively as follows: Expanding 



(4.6a) and (4.6b) and assuming for simplicity that a{k) has no zeros, we find 
2 



9u 



111 



kxn{t,k) - ]dk 



2 

ni 
4 



[{a^ — a)bi{ak) + (a — a^)6i(a^A;)] kt 



dk, 



521 



+ 



4 



k'^Xii{t,k) 



3fcgoi \ 
2 J 



dk 



kh-^'^''[{l - a)bi{ak) + (1 - a'^)bi{a^k)]dk. 



(4.18a) 



(4.18b) 



Using equation (4.17a) to determine xiij we can determine ^n, g2i from (4.18); then 
^11 can be found from (4.16). These arguments, just like in the case of NLS, can be 
extended to higher orders and thus yields a constructive scheme for computing A{k) 
and B{k) to all orders. 

Similarly, for the first Neumann problem, the expansion of (14. 71) gives 



9oi{t) 

92l{t) 



- [ Xndk - - 

JdDs ^ JdD3 



[(a — a^)bi{ak) + (a^ — a)bi{a^k)\dk, 



lXu{t,l)dl dk 



(4.19) 



4 

+ - 

TT 



^2g-8ifc3t 



[(1 - a^)bi{ak) + (1 - a)bi{a'^k)]dk. 



In this case, we first use equation (4.17b) to determine xii; we then determine 511, 521 
from (4.19); then <I>ii can be found from (4.16). For the second Neumann problem, 
we use a similar argument based on (4.17c). In all cases, the system can be solved 
perturbatively to all orders. 

4.4. The mKdVII equation. An important question in the analysis of initial- 
boundary value problems is the determination of the number of boundary conditions 
needed for a well-posed problem. For the mKdV equation (2.4a) considered earlier 
(in this subsection we will denote this equation by mKdVI), we have seen that one 
boundary condition must be imposed at t = 0. On the other hand, for the mKdVII 
equation, which has a minus sign in front of d^q/dx^, 



dt dx^ dx 



0, 



q real, A = ±1, 



(mKdVII) 
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two boundary conditions must be prescribed at t = for a well-posed problem. 
How does this difference between mKdVI and mKdVII manifest itself in the present 
formalism? 

To answer this question we will assume for simplicity that qq = 0. There exist 
three different representations for go{t) for mKdVI: 

9o{t) =^ [ Xidk, goit) =^ [ xidk, go{t) = - [ xidk. (4.20) 

The latter two representations are useful for the effective solutions of the first and sec- 
ond Neumann problems, respectively. Indeed, the perturbative equations (cf. equa- 
tion ( |4.17| )) 



rt 

Xinit, k) = 12/c2 / e'^*'(*'^*)5on(t')'^*' + lower order terms, (4.21a) 

rt 

X^^[t, k) = 6ik / e^*^'(*'"*)5i„(t')rfi' + lower order terms, (4.21b) 
Jo 
j-t 

X^^{t, k) = -3 / e^'^'^*'~^^g2n{t')dt' + lower order terms, (4.21c) 
Jo 

show that if gi and g2 are known, then xi and xi respectively can be determined 
uniquely at each step of the perturbative expansion. 

The second representation in (4.20) is derived by noting that 



f Xidk= [ +a W / )^idk. (4.22a) 

Since is bounded and analytic in C \ D^, each of the integrals on the right-hand 
side can be computed individually and this yields 

Xidk = — a— - a — ^ = irgo- 



dDs \ 3 3 3 , 

The third representation is derived in the same way by noting that 

f x^dk=([ Wf +a [ ]<^>idk. (4.22b) 

Let us now consider the situation for mKdVII. We claim that the representations 
in (4.20) have the following analogs for mKdVII: 

9oit) = -- [ Xidk, go{t) = -- [ {xi+Xi)dk. (4.23) 

^ JdD3 ^ JdD-i 

Indeed, mKdVII admits the Lax pair 

f//^ - iA;[cr3,/i] = g/x, ^^^^ 



where 



Q(x, t) = ( ^"j , g(x, t, k) = -2Q3 + Q^^ _ 2ik{Q^ + Q,)a3 - Ak'^Q. 



Proceeding as in the case of mKdVI, we find 

2ik ^ U2. 
this relations] 

in the case of mKdVI, the first representation in (|4.23l). Morever, equations (|4.22l) 



i.e. go = —2i^^\ Using this relationship it is straightforward to derive, just like 
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are also valid for mKdVII. However, in this case the analyticity domain C \ of 
<I>i is smaller (see figure 6]) and this prohibits the evaluation of the integrals on the 
right-hand sides of (4.22). In order to proceed, we now have to consider the sum of 



equations (4.22a) and (4.22b): 



9^3 



(Xi +Xi)dk 



(1) 
1 ' 



which yields the second representation in (4.23) 



In summary, whereas there exist three different representations (4.20) of go for 



mKdVI, there exist only two such representations for mKdVII. Since the perturbative 



equations (4.21) are also valid for mKdVII provided that the signs of gon and g2n are 



reversed, the second representation in (4.23) can be used to determine go at each step 



of the perturbative scheme provided that both gi and 52 are prescribed. 

Similar remarks apply to the representations of gi and 52- For example, for (72 
there exist three different representation for mKdVI involving xi, xi and xi! respec- 
tively; the former two are useful for the effective solutions of the Dirichlet and first 
Neumann problems, respectively. On the other hand, for mKdVII there exist only 
two representations involving xi + Xi and xi respectively (cf. equation (1.18a) in 
|35j): the former can be used to give an effective solution provided that both go and 
gi are prescribed. 

5. The Dirichlet to Neumann map 

We have shown in sections [3] and |4] that the spectral functions associated with the 
Dirichlet and Neumann problems for NLS and mKdV are characterized by nonlinear 
integral equations which can be solved perturbatively to all orders. In this section, 
we demonstrate that these integral equations can also be employed to construct per- 
turbatively the generalized Dirichlet to Neumann map. 

Consider the Dirichlet problem for NLS on the half-line. The perturbative con- 
struction of the Dirichlet to Neumann map involves the following: Given initial data 
qo{x) and Dirichlet data go{t) in the form 

qo{x) = eqoi{x) + e^go2(x) H , goit) = egm{t) + e^gQ2{t) H , 

where e > is a small parameter, determine the coefficients {gin{t)}n>i in the ex- 
pansion of the corresponding Neumann data, 

gi{t) = egii{t) + e^gi2{t) + e^g^sit) + • • • , (5.1) 

in terms of {qonit), gon{t)}n>i- Assuming for simplicity that qo{x) = 0, we will derive 
expressions for ^n, 512, and gi3 in terms of goi, go2i 903 for NLS on the half-line. 
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Similar results can be obtained for the Neumann to Dirichlet map for NLS, as well 
as for the various generalized Dirichlet to Neumann maps for mKdV. 

Theorem 5.1. Let 

q{x, t) = eqi{x, t) + 6^^2(2;, t) -\ , 

be the perturbative solution of the NLS equation on the half-line satisfying the initial 
conditions q{x, 0) = 0, x > 0, and the Dirichlet boundary conditions 

q{0,t) =egoi{t) + e^go2{t) + -- - , 

where {gonit)}n>i o-fc sufficiently smooth functions compatible with the zero initial 
data, i.e. goniO) = ffon(O) = ffon(O) = • • • = 0, n > 1. Then, the first few coefficients 
in the expansion (5.1) of the Neumann data gi{t) = qx{0,t) are given by the following 
formulas: 



guit) 

513 (i) 



e 4 



^^dt , g^t) 
2Ac f'\9oi{t')\^9oi{t')^. 2Ac 



e 4 



Vt - t 



TTl 



Vt^ 

2c r 503(t') 



dt' 



\9oiitT 



+ 



+ 



vr 

Ac 

iri 

Ac 

IT 

Ac 
vr 

Ac 
vr 



501 (i" 



(5.2) 



Vt^ 



goi{t') / 5'ii(*") 



511(^0 / goiit") 







t' 







gu{t') I -gu{t") 
■t 

'0 



{t-t' + 1" - 1'")^/^ 

gain 
{t-t' + 1" - 1'")^/^ 

gain 
{t-f + 1" - f'fl'^ 

golin 



{t-f + 1" - f'fl'^ 

+ iXgoiit) I [goi{t')giiit') - gii{t')goi{t')]dt', 



dt"'dt"dt' 
dt"'dt"dt' 
dt"'dt"dt' 
dt"'dt"dt' 



(5.3) 



y. — ^ A 

2 



Proof. Let 



^>i = e$ii + e^$i2 + 



$2 = e^>2i + e^^>22 + 



and let {Xjli ^"^^ {Xjli be the odd and even combinations formed from {•J'jli as in 
(|3l2| and $1^. Let 



Xj = eXji + e Xj2 + 



2 - 



We first prove (|5^. By (|3.30|), 

9ii{t) 



2 

Tii 



dD3 



Xj = eXji + e Xj2-\ 



{kxn{t, k) + igoi{t))dk. 



(5.4) 



Substituting the expression (3.31) for xii into (5.4) and integrating by parts, we find 

rt 



guit) = - , 

JdDi Jo 



goi{t')dt'dk. 



(5.5) 
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Using the identity 



vr 



yt 2 



(5.6) 



we can compute the k integral in (5.5). This gives the first equation in (5.2). The 
derivation of the second equation in ( 5.2[) is s imilar. 

In order to prove (5.3), we note that ( 3.11| ) implies that 

2 



513 (t) 



{kxUt,k)+igomdk+'^-^^ 



9^3 



X22{t,k)dk. (5.7) 



Moreover, by (|3.10) 

$13 



f e^ikHt'-t) [-iAIsoiP^-ii + 2A;go3 + ^513 + (2%oi + %i)«'22] dt', 
Jo 

$22 = A / [(2A;5oi - %i)$ii + iboiP] dt' , 
Jo 



so that 



Xl3 



AikHt'-t) [_iA|5oiPxii + 4%o3 + 2%oiX22 + ^511X22] dt' , 

X22 = A / [2A;5oiXii - %iXii] dt' , 
Jo 

X22 = A / [2/0501X11 - %iXii + 2i|5oin dt', 
Jo 



Xii 



2i I e^"'"^''-'^gndt'. 



I 



Thus, letting / denote the first integral on the right-hand side of (5.7), i.e. 

{kxri{t,k) +igos{t))dk, 



we find 
I -- 



/ e 

dD3 I Jo 



4jfe2(t'-t) 







-iX\goi{t')\Mk I e 
ft' 

'0 



ft' 

+ 4kgo3{t') + 2kgoi{t')\ / [2/cgoiXii - %iXii + 2^1501^] (i", k)dt" 

Jo 



mi (*') A / [2A;goiXii - %iXii] (*", k)dt" 







dt' + igo3{t) >dk 



^A / / |5oi(i')P / 4k^e^'^'^'('"-''>goi{t")dt"dt'dk 



+ 



+ 



dD-i Jo 



^.^"'"^''~'hk^go3{t')dt' + igo3{t)] dk 



OD3 Jo 



,^^k-it'-t)2k'g^^(t')\ 



goi{t"W 



"^8i.2 / .^^k'ii"'-ng^^^t"')dt"' 



i-gu{t")2i / e'^'''^'"'-'"'>gii{t"')dt"' + 2i\goi{t")\' 



dt"dt'dk 
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IdDa Jo 



2kgoi{t")2i / e 



^4ik'^(t"'-t") (.III 



9ii{r)dt" 



dfdt'dk. 



Integration by parts yields 

/ /' 

IdDi Jo 



iX 



\9oi{t')\H-ie'''"'^''-'^goi{t')+i I e'^^'^'"-'^m{t")dt"yt'dk 



+ z e 

IdD-i Jo 



dD3 



-igoi{t)X 



dosit )dt dk 



-goi{t')i / e''''^'"'''^goi{t")dt" 



dt' 



+ -giiit") [ 

Jo 



3'^'''^'"'-'"^gn{t"')dt"' 



dt" Ut' Uk 



[ f e'^^'^''-%g^,{t')X f Li(i") f - f " e^^^^(*"'-*'\n(t'")dt" 

JdDz Jo JO L V Jo 



gii{t"){ goi{t") 



,^^kHt"'~t")g^^(^tiii^dtiii 







dfdt'dk. 



Using the identity (|5.6p to perform the dk integrals, we deduce that 

goijt" 

'0 \V^ — ''^' Jo Vt — t-- J Jo 



I = Xc \goi\ 



;( goijt') 
\Vt-t' Jo Vt - 1" 



mit') 



+ iXcgQi{t) 



iXc / goiit') 



igoiit') 
t' 





t' A f+ll 



golit" 



-M" + guit') 







igoi{t") 



501 (t" 



Vt-t 



df 



gii{t" 



Vt-t' + 1" - 1'" 



Vt^^ 
dt'" 



dt" 



dt' 



Vt-t' + t" - t' 



.dt" 



dt"dt' 



-iXc / goi{t') 



igoiit') 



t' X (+11 



goi{t" 



,dt" + gu{t') 



t' „ r+ii 



t „ ui\ rt' 



tXc I f [m{t'')gii{t'')--g,,{t'')goi{t'')\dt"dt' 

Jo yt — t Jo 

+ f gu(t') f gm(t") I" 



sun,,, 

Vt^ 



dt' 



--dt"'dt"dt' 



iXc 



fguit') f 

Jo Jo 



-gu{t") 



Vt-t' + t" - t' 

goijt'") 
Vt-t' + 1" - 1' 



-.dt'" dt" dt' . 
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We next integrate the terms involving 1/ ^Jt — t' + 1" — t'" by parts with respect to 
dt' . The total contribution of lines 2-5 of the above expression yields 



iXc / gm{t') 



mi{t") 

+ -9ii{t" 



goiit" 



2{t -t' + t" - t"'f/^ 

giiit'") 



dt" 



-dt" 



dt" dt' , 



2{t -t' + t" - t«<fl'^ 
while the last two lines combine nicely with the sixth line. We thus arrive at 



I = Xc \goi 



9oi{t') f mit") 



yjt - t' Jo Vt - t" 



zdt" \ dt — ic 



Ac / goi{t') / -goi{t") 



mit" 



+ iXc / goiit') / giiit") 



Ac / ign{t') / goi{t") 



+ Ac / igii{t') / gii{t") 
Jo Jo 

On the other hand, using the identity 



k / e^'''"^'"-''^f{t")dt"dk 

dDz Jo 



2{t- 


- t' + t" - 


t'")3/2 




9u{t"' 




2{t 


-t' + t" - 


- t'")3/2 




9n{t"' 




2{t 


-t' + t" - 


- t'")3/2 




goi(t'" 




2{t 


-t' + 1" - 


_ ^///-)3/2 



'^dt' 

Vt^ 
it"'dt"dt' 

dt'"dt"dt' 
dt'"dt"dt' 
dt"'dt"dt'. 



(5i 



f/(t'), 0<t'<t, 



0<t' = t, 

where f{t") is a smooth function, we find that the second integral on the right-hand 
side of (5.7) can be written as follows: 

r-t r rt' 



X22{t,k)dk 



X 

dDz Jo 



2kgoi{t')2i / e'"^'^'"-'">gu{t")dt" 



igu{t')Ak [ e 
Jo 



^ 'goi[t )dt 



dt'dk 



iX- / [goi{t')9ii{t')-9ii{t')goi{t')\dt'. 
^ Jo 



(5.C 



Equations (5.7), (5.8), and (5.9) imply (^5.3|. 



□ 



6. Conclusions 

As it was mentioned in the introduction, linear PDEs involving only second order 
derivatives can be analyzed by an appropriate extension from the half-line to the 
full line. In analogy with the linear case, linearizable BVPs for the NLS and the 
sG equations have been studied via techniques based on an appropriate extension 
from the half-line to the real line. These extensions yield explicit conditions on the 
scattering data for an initial value problem formulated on the full line, see [U [5l [28l 
[51182]. 

It is of course possible to map the RH problem obtained by the unified method to 
the RH problem obtained via an extension from the half-line to the full line, see for 
example [36]. However, as it was noted in the introduction, the unified method has 
the advantage that it always yields solutions that are uniformly convergent at the 
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boundaries. Furthermore, the unified method is apparently the only approach that 
can be apphed to hnearizable BVPs for PDEs involving third order derivatives such 
as the KdV and the mKdV equations. 



Appendix A. A comparison of the unified method with the 

APPROACH of FTHI 



In sections |3] and |4] an effective characterization of the generalized Dirichlet to 
Neumann map for the NLS and the mKdV equations respectively, was presented. 
The effectiveness was demonstrated by showing that the relevant nonlinear equations 
can be solved perturbatively via a well-defined recursive scheme. Furthermore, in the 
linear limit, the relevant formulas coincide with the formulas obtained by solving, via 
the new method, the associated linearized equation. 

It was noted in the introduction that for the special case of the NLS on the half-line 
(but not for the mKdV on the half- line, or even for the NLS on the finite interval), 
the approach introduced in [18j yields expressions for go{t) in terms of qo{x) and gi{t) 
for the Neumann problem and similarly for the Dirichlet and the Robin problems. 
In what follows we compare the unified method with the approach of [18] for the 
following slight generalization of the linearized version of the NLS: 



iut{x, t) + u^xix, t) + iaux{x, t) = 0, 
where a is a real constant. 



0<a;<oo, 0<t<r, 



(A.l) 



A.l. The unified method. We first implement the new method. Equation (A.l) 
can be written in the following divergence form: 



ue 



-ikx+iw{k)t 



0, 



k G 



where the dispersion relation w{k) is given by 

w{k) = + ak. 



(A.2a) 



(A.2b) 



We note that equation (A.2a) and Green's theorem immediately imply the global 
relation 



J^ik)tu(^k,t) = UQ{k) -~g{k,t), Imk < 0, 



(A.3) 



where ■u(/c, t) denotes the half-Fourier transform of u{x, t) and the spectral functions 
{uQ{k),g{k,t)} are defined as follows: 



UQ{k) 



uo{x)e ^^"^dx, 



k £ 



(A.4) 



g{k, t) = igi{w{k),t) - (a + k)gQ{w{k),t), k £ C, 



9j{k,t) 



j = 0, 1. 



(A.5a) 
(A.5b) 



We note that equation (A.3) can also be obtained by the usual application of the 
half-Fourier transform. 



Employing the inverse Fourier transform, equation (A.3) yields 
1 



u{x, t) 



2ti 



-ikx—iw{k)t f 



UQ{k)dk 



1 

2^ 



-ikx—iw{k)t ~ 



g{k,t)dk. (A.6) 
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Jordan's lemma implies that it is possible to replace in the second term of (A. 6) the 
integral along the real axis with an integral along the contour dD which consists of 
the union of two rays: 



dD : 



Imk = 0,-^ <Refc<ooJ>U<!Re/c 



a 



, < Im k < oo 



We next consider the transformations k — t- i^(A;), which leave w{k) invariant: 
k^ + ak = u'^ + au, i.e. u = k and u = —k — a. 
Replacing in equation ( A.3| ) k with v{k) = —k — a, we find 



Jw{k)t 



u{—k — a,t) = uo{—k — a) — igi — kgo = 0, 



Im > 0. 



(A.7) 



In the case of the Neumann problem, we use (A.7) to eliminate i-e- to eliminate 
the transform of the unknown Dirichlet boundary values. Solving equation (A.7) 
for go, substituting the resulting equation in (A. 6), where the contour in the second 
integral is given by dD, and employing Jordan's lemma to show that the contribution 
of u{—k — a, t) vanishes, we find 



u{x, t) 



1 

2^ 



-ikx—iw{k)t 



uo{k)dk 



1 

2^ 



-ikx—iw(k)t 



dD 



i[2+-]gi{w{k),t) 



a 



1 + ^ ]uo{-k-t) 



(A. 
dk. 



Remark A.l. 1. Equation (A. 8) expresses u{x,t) in terms of the Fourier transform 
uo{k) of the initial data uo{x) and of a t-transform of the Neumann data gi{t) := 
dxu{0,t). 

2. By employing Jordan's lemma, it follows that gi{w{k),t) can be replaced with 
gi{w{k)), where 



Jw{k)T 



gi{T)dT, k £ C 



A. 2. The approach of [18]. We next consider the approach of 
case of a = 0. In this case, equation (A.6) becomes 

u{x, t) 



1 

27r 



-ikx—iw(k)t 



1 

2^ 



UQ{k)dk 



[i~gi{w{k),t) - kUw{k),t)] dk. 



(A.9) 
in the special 

(A.IO) 
(A.ll) 



The important new idea introduced in [181 has the following implementation for linear 
PDE s: Eli minate directly go from equation (A.IO). In this respect, evaluating equa- 
tion (A.IO) at = 0, multiplying the resulting equation by exp(iA:^T), and integrating 
with respect to r from r = to r = t, we find 



1 

2^ 



Jk^ 



^uo{l)dl 



^l~gi{l\T)dl-h{T)\dT, 



where 



h{r) = / 



le 



gQ{f,T)dl. 



(A.12) 
(A.13) 
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Due to the invariance of go{P,T) under the transformation I 
and (A. 12) becomes 

1 

2^ 



I, loiT) vanishes 







uo{l)dl 



igi{f,T)dl} dr. (A.14) 



This equation express es dir ectly gp in terms of the known transforms uq{1) and 
gi{P, t). Substituting (A.14) into (A. 10) we find a formula expressing u{x, t) in terms 
of the given initial and boundary conditions (albeit, this formula is more complicated 
than equation (A. 8) evaluated at a = 0). 

It appears that the above approach cannot be applied to PDEs whose dispersion 
relation is not invariant under the transformation k ^ —k. For example, if a 7^ 0, 
the analogue of equation (A. 12) is now the following equation: 



)iwik),t) 



1 
2^ 







-iw{l)T 



uo{l)dl 



T)d/ - I(r) \dT, (A.15) 



where 



I(r) 



-'-'(^)-{a + irgo{w{l),r)dl. 



In this case I(r) 7^ 0, hence (A.15) does not provide an explicit expression for ^0 in 
terms of the given initial and boundary conditions. 

The above discussion suggests that the effective treatment of linear and integrable 
nonlinear PDEs requires the explicit use of the transformations which leave the lin- 
earized dispersion relation invariant. 

Appendix B. The asymptotics of c{t,k) 



3.1 



the proof of lemma 4.1 is similar. The functions {^jlf 

(B.l) 



We will prove lemma 
satisfy 

= -Aik'^^i - iAIffoP^i + (2%o + igi)^2 
$2* = A(2A;5o - igi)^i + M\go\'^'^2- 
It follows that these functions admit an expansion of the form (see chapter 6 of [1 



^i{t,k) 
'^2{t,k) 



ao(t) + "# + 



+ (/3o(t) + ^ + 



k 



k 



e 

00, k £ 



(B.2) 



where the coefficients aj(t), /3j{t), j > 0, are column vectors which are independent 
of k. We determine the coefficients by substituting (B.2) into (B.l) and using the 
initial conditions 



ao(0) + /3o(0) = (0, 1)^, ai(0) + /3i(0) = (0, 0) 



T 



This yields 

Mt,k) 
^2{t,k) 



+ 



+ 



$^(0) 



+ 



.(2) 



k- 



00, k £ 



(B.3) 
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Substituting these expansions into the global relation (3.2), we find 

,,,, . £l!w . ££V2 . of J,) . f ^iii^ . of J,) V--. (B.) 



The assumption that c{t,k) is of 0{l/k) as /c — )• oo in C"*" implies that the terms in 
(B.4) involving e"^*'^ * must vanish, i.e. for consistency we require 

6(1) = $S^)(0). (B.5) 



Using the expressions for b^^^ and $|^^^(0), we find that equation (B.5) is valid ifi' the 
initial and boundary conditions are compatible, i.e. iff 

go(0) = 50 (0). 

Example B.l. When A = —1, NLS admits the one-soliton solution 

q{x,t) = e*''sech(x - 2t). 
For this solution it is possible to compute {^j}i explicitly: 

e-^*'^'* ( i{2k + 1) f _2e2*+4^'='* + e^* + l) - e^* + l' 
<^>iit,k) = ^ 



2g2t-4ifc^t 
<l>2{t,k)=^^' 



((2A; + 1)2 + 1) (e4t + 1) 

^^t^ + {2k + l){2k - itanh(2t) + 1) 



(2A: + 1)2 + 1 

It is easy to verify that the relevant asymptotics is of the form ( |B.3 ) . 

Appendix C. The necessity of computing A{k) and B{k) 

We consider NLS with vanishing initial conditions and 'small' Dirichlet boundary 
conditions, 

g(0,t) = egoi(t)+e'5o2(t) + e'go3W + 0(e'), e ^ 0. (C.l) 
We expand the unknown Neumann boundary values in the form 

g,(0,t) = e(7ii(t) + 6^2(0 + eV3(t) + 0(e'), e ^ 0. (C.2) 



The global relation (2.17) becomes 

A{k) = F2{T,k), B{k) = -e^^*='^Fi(r, k), k G C+, (C.3) 

where Fi and F2 are analytic functions of k in C+ and of order 1 + 0(1/A;) and 0(1 /k) 
respectively as /c —t- 00 in C^. 

The spectral functions can be expanded in the form 

A{k) = l + e^A2ik) + 0{e^), 

B{k) = eBi{k) + e^B2{k) + e^B:i{k) + O(e^), e ^ 0, /fc G C. (C.4) 
The definition of B{k) yields 

Bi{k) = i r e^'^^^ gii{T)dT - 2k f e"^'^'^ goi{T)dT, kGC. (C.5) 
Jo Jo 

This implies the important symmetry relation 

Bi{-k) = Bi{k) + Ak e^''''^goi{T)dT, k£C. (C.6) 
Jo 



34 



A. S. FOKAS AND J. LENELLS 



Replacing in the 0(e) terms of the global relation (C.Sa) k with —k and then using 
the symmetry relation ( |C.6 ), we find 



Biik) 



-Ak 



T 



'goiiT)dT 



>i(r,-A;), ke 



(C.7) 



The function Fi(T,—k) is unknown, however using the fact that it is bounded and 
analytic in k for k G C~, it follows that its contribution vanishes. 

It turns out that the situation for B2{k) is similar, however, the unknown functions 
entering in the right-hand side of the equation for B^[k) yield a nontrivial contribu- 
tion. Indeed, in order to express {A{—k), B{—k)} in terms of {A{k), B(k)} and go{t), 
we replace k by —k in the ODE (2.13) satisfied by denoting by " the operation of 
replacing k with —k we find: 



$t + 2ik^a3^ 



-2kQo + Q''^^)^. 



Hence, equations (2.13) and (C.8) imply 



4ke 



The (12) component of this equation yields 
^BA - AB 



Ak ( e'^'^^^goAA - Xe-^'^''^goBB 



k £ 



(C.8) 



(C.9) 



(C.IO) 



The 0(e) term of this equation yields equation (C.6), whereas the O(e^) and O(e^) 
terms yield the following equations: 



B2 = B2 + 4:k 



(C.ll) 



S3 = + Ak 



go3{T)dT + U{k), keC, 



where 



T ^ 



U{k) = B1A2 - B1A2 + 4fc / e'^''' "501(^2 + A2)-e 



-4jfc^T 



goiBiBi 



(C.12a) 



dr. (C.12b) 



Replacing in the O(e^) and O(e^) terms of the global relation (0.3^) k with —k and 
using equations ( |C.ll ) and ( C.12| ), we find 



B2 = -4k I e"'^ ^gm{T)dT - e^'^ ' F2{T, -k) 
3 



/fceC", (C.13) 
B3 = -4kl e''^^^go3{T)dT-e''''^' F3{T,-k)-U{k), k e £- . (C.14) 







The contribution of the terms involving F2{T, —k) and F^iT, —k) vanishes, however, 
U{k) involves terms, such as which are not bounded for k G C~, thus the contri- 
bution from U{k) does not vanish. This makes it necessary to determine Bi in terms 
of 50- 
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Appendix D. A less effective approach 



It is possible to characterize Fi{t, k) and F2{t, k) in terms of a system of nonlinear 
singular integrodifferential equation. Indeed, the functions F and F satisfy equation 
(C.9) but only for k gM. The (12) component of this equation yields 



e4ifc2i ( ^^^^ _ p^p^ 



Ake 



k £ 



^doFiFi — \goF2F2 

Assuming that Fi{t, k) has no zeros for k € C^, the determinant condition 

F2F2 - XFiFi = 1, keR, 

yields 



(D.l) 



F2{t,k) 



2^/r^ln(l+|FlP(t,/)) 



dl 
l-k 



k£C-^ 



(D.2) 



The function Fi{t, k) is bounded and analytic for A; G C , hence it can be represented 
in the form 



Fi = f + iHf, Fi = f- iHf, Fi = f- iHf, 



where 



iHf){k) 



nT_r^ I - k ' 



k £ 



(D.3) 



(D.4) 



Substituting equations ( |D.2 ) and (D.3) in equation (D.l) and equating the real and 
imaginary parts of the resulting equations, we obtain two equations for the two un- 
known functions / and /. 



Remark D.l. 1. The (11) component of equation (C.9) yields 



FiFi 



XF2F2 )^ = 4AA; ( goFiF2 - goFiF2 



k G 



Each term of this equation is bounded and analytic for k G C , thus the substitutions 



(D.3) yield a single equation, which actually turns out to be an identity. 



2. The appearance of the Hilbert transform makes the above system of equations 
hard to analyze. A system of singular integrodifferential equations was also presented 
in yj], see equation (4.6). Actually the first such system of equations was presented 
in [28] by using an odd extension from the half- line to the full line. It appears that 
the system of [28] is simpler than both the system presented in [19] and the system 
presented here. 
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